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For controlling periodic orbits with delayed feedback meth¬ 
ods the periodicity has to be known a priori. We propose a 
simple scheme, how to detect the period of orbits from prop¬ 
erties of the control signal, at least if a periodic but nonva¬ 
nishing signal is observed. We analytically derive a simple 
expression relating the delay, the control amplitude, and the 
unknown period. Thus, the latter can be computed from ex¬ 
perimentally accessible quantities. Our findings are confirmed 
by numerical simulations and electronic circuit experiments. 

PACS numbers: 05.45.-|-b, 02.30.Ks, 07.50.Ek 

Control techniques using time-delayed output signals 
are a very well established field and known for at least 
half a century in the engineering and mathematical con¬ 
text (e.g. 0 and references therein). Delayed feedback 
control methods, which have for the physicists’ purpose 
been rediscovered in , are very useful since neither spe¬ 
cial knowledge of the system under consideration nor so¬ 
phisticated reconstruction techniques are required, and 
the method is easily implemented in experiments . As 
a certain kind of drawback, the success of delayed feed¬ 
back methods is difficult to predict, and the stability 
analysis of the corresponding delay systems shows a rich 
behaviour (e.g. Q). Only recently some progress in the 
understanding of general features has been made in the 
physical context Q. Since control of actual periodic or¬ 
bits with delayed feedback methods requires a delay time 
which is an integer multiple of the period, one runs into 
principle difficulties whenever the period is not known a 
priori. Some empirical schemes have been reported to 
circumvent such problems |^. They work quite well for 
special cases but no theoretical foundation has been pro¬ 
posed. Here we address the problem that the period of 
the unstable periodic orbit is unknown. A systematic 
strategy is developed to obtain the desired period, when¬ 
ever a periodic control signal is observed. 

Theoretical approach - To keep our approach as general 
as possible the theoretical considerations are based on a 
fairly arbitrary equation of motion 
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X = F{x{t),K{g[x{t)]-g[x{t-T)])) . (1) 

Here x denotes the phase space variables, g[x] the mea¬ 
sured scalar quantity, r the delay time, and K the control 
amplitude. We do not specify the functional dependence 
of the systems on the control signal g\x{t)] — g[x{t — r)], 
since this dependence is in general difficult to estimate 
from the experimental point of view. Without control, 
AT = 0, the system should admit an unstable periodic 
orbit ^it) with period T and Floquet exponent A -|- iw, 
A > 0. We intend to stabilise this orbit. 

Whenever the delay differs from the period, t ^ T, 
the orbit ^ does not yield a solution of the system sub¬ 
jected to control. However, the system admits a periodic 
solution T] with period 0. Such a statement can even be 
proven rigorously 0 provided that the delay mismatch 
r — T is not too large. In addition, the fictitious solution 
r] tends towards the unstable orbit ^ in the limit t ^ T. 
Of course, the period of this fictitious orbit depends on 
the parameters of the system, in particular on the de¬ 
lay time and the control amplitude, 0 = Q{K,t). We 
remind the reader, that the quantity 0 can be observed 
from the period of the control signal, whenever the orbit 
T] is stable. In what follows we assume that the sys¬ 
tem parameters are adjusted in such a way, i. e. we can 
observe the period 0 for different values of the control 
amplitude K and the delay time r. 

The strategy for the determination of the desired pe¬ 
riod T is quite simple. Since the orbit ^ yields a peri¬ 
odic orbit of the controlled system for r = T, the mea¬ 
sured period of the control signal obeys Q{K,T) = T. 
Hence we simply have to look for zeros of the function 
Q{K,t) — T. The latter can be measured in principle, 
provided we meet the assumption made above. Neverthe¬ 
less, it would be helpful if some analytical result about 
the dependence of 0 on the delay and the control am¬ 
plitude would be available. We show that up to second 
order in the mismatch t — T the relation 

0(iT,r)=r+-^(r-T) + O((r-T)2) (2) 

holds. Here k denotes a system parameter which cap¬ 
tures all the details concerning the coupling of the con¬ 
trol force to the system. Since the parameters r and 
K are adjustable in experiments and 0 is a measurable 
quantity, the desired period can be computed from eq. (|^) 
using two data points. 
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In order to derive expression we rewrite eq. (0) for 
the periodic orbit rj in terms of the dimensionless time 
s = t/Q as 

r]'{s) = eF{q{s),K{g[r]{s)] - g[r]{s - r/0)])) (3) 

and 

vis) = vis+ 1) ■ (4) 

Since eq. (i represents a Taylor expansion we are looking 
for the derivative cli-0|r=T- For that reason one takes the 
derivative of eq. i) with respect to r, keeping in mind 
that the periodic solution rj depends explicitly on r. 

idrV)' - QDiF{. . .)drVis) - &Kd2F{. . .) 

■ {Dg[vis)]drVis) - Dg[Tj{s - T/e)]drTj{s - t/Q)] 

= idrQ)F{...) + QKd2F{...) 

■ {Dg[Tiis-T/Q)]-n'{s-T/Q)}driT/Q) . (5) 

Here Di and d 2 denote the derivative with respect to the 
first/second argument of F , and the arguments abbrevi¬ 
ated by ■ • ■ coincide with those from eq.(^). The contri¬ 
butions involving the derivative of the orbit with respect 
to the explicit r-dependence, drV^ have been collected 
on the left hand side. The boundary value problem (||), 
(P determines both, drQ as well as drV- In order to 
separate the former quantity we trace back to the fact 
that the linear operator on the left hand side of eq.(H) 
admits a vanishing eigenvalue. The corresponding Gold- 
stone mode is related to the translation invariance in time 
of the original system. In fact, taking the derivative of 
eq.(^) with respect to s one obtains 

0 = (rj')' - QD^F{.. .)ri'{s) - QKd 2 F {...) 

• {Dg[vis)]v'is) - Dg[r]{s-T/Q)]r]'{s-T/Q)} . (6) 

Eq.(^ just states that t]' yields the right-nullei- 
genfunction. Within the canonical scalar product 
/p v{s)u{s)ds we denote the corresponding left-nullei- 
genfunction by C('S)- All the terms on the left hand side 
of eq.(||), which involve OtT] vanish identically after mul¬ 
tiplication with (■. Hence we are left with 

O = a^0 f Cis)Fi...)ds + QKdriT/Q) 

Jo 

■ ( Cis)d 2 F{...){Dg[q{s-T/Q)]'q'{s-TlQ)}dail) 
Jo 

The details of the system, which are only contained in 
the integrals, are now condensed to simple numbers. But 
in general the integrals depend on the delay r and in 
particular on the control amplitude K through the left- 
eigenfunction C (cf. eq.(^)). For that reason we evaluate 
eq.(0) at T = T. Then 0 = r holds and the delay in the 
arguments of rj drops by virtue of the boundary condition 


(^. Due to the same argument the linear operator (||) 
and therefore the eigenfunction becomes independent 
of K. Hence the integrals become constant real numbers 
and eq. (0) yields 

0 = KdrQ\r=T+TKdriT/Q)\r=T (8) 

Here n denotes the ratio of the integrals occurring in 
eq. ( 0 )- We solve for dTQ\T=T^ and obtain eq.(^ from a 
simple Taylor series expansion. 

Numerical simulations - We demonstrate the applica¬ 
bility of our analytical results by numerical simulations 
in an autonomous system. First of all stabilisation of 
periodic orbits by delay methods requires a finite tor¬ 
sion, i. e. a finite frequency in the Floquet exponent of 
the controlled orbit |^. Since autonomous equations al¬ 
ways admit a vanishing exponent a finite frequency can 
be realised in dissipative three dimensional models only 
by a complete flip of the neighbourhood of the orbit. For 
that reason certain equations like the Lorenz model can¬ 
not be stabilised at all by delay methods, apart from the 
fixed points for which the reasoning given above does 
not apply. Therefore we concentrate here on the Rossler 
equations as a certain minimal model for our purpose. 

ii = -X 2 -x^-K ig[x{t)] - g[x{t - t)]) 

X 2 =xiF ax 2 - K ig[x{t)] - g[x{t - r)]) 

xz = h + xixz - cxz (9) 

Our results do not seem to depend significantly on the 
coupling of the control force to the original equations of 
motion and on the particular choice of the scalar quan¬ 
tity g{x]. We have used a bounded quantity in order 
to avoid diverging solutions. The results presented here 
correspond to the choice g[x\ = tanh[(a;i -I- a;2)/10]. In 
addition, the system parameters have been fixed to the 
values a = 6 = 0.2,c = 5.7to ensure chaotic dynam¬ 
ics in the absence of control. For our control purpose 
we concentrate on the period-two orbit in the canonical 
Poincare map with T = 11.758 ..., XT = 1.256 ..., and 
ujT = TT. Numerical simulations have been performed 
by means of an adaptive stepsize Runge-Kutta method 
of order four, together with a cubic spline for the delay 
from the NRF library [||. 

For a quite large range of delay times r one observes 
two critical values of the control amplitude which limit an 
interval where a stable periodic orbit r] can be observed. 
From the Fourier transform of the scalar quantity g[x\ it 
is evident (cf. fig.0) that at the lower critical value the 
orbit loses stability via a flip bifurcation, whereas at the 
upper critical value a Hopf bifurcation occurs. In order 
to check the accuracy of eq.(||), the period 0 of the fic¬ 
titious orbit has been extracted from the peaks in the 
Fourier spectra of the control signal. The dependence of 
0 on the control amplitude for several delay times is sum¬ 
marised in fig]^ and compared with our analytical expres¬ 
sion. The apparent systematic deviation of the analytical 
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result just comes from the fact that the latter is a hrst- 
order approximation to the curved manifold Q{K,T) in 
the three-dimensional K-T-Q space. In summary, eq. (|^) 
describes the observed periods quite accurately. 




FIG. 1. Absolute square value of the Fourier transform of 
the scalar quantity g[x\ for r — T = 1.0 in the vicinity of the 
lower and upper stability threshold. The transform has been 
performed for a series of length 1024 x r discarding a transient 
of 100 X r. The spectrum has not been properly normalised 
and the abscissa extends over 9 decades. 


Experiments - To illustrate the experimental accessi¬ 
bility of our analytical results we have performed mea¬ 
surements on a nonlinear electronic circuit (cf. hgj^). The 
circuit consists of several operational amplifiers (three 
acting as integrators, two as inverters) with associated 
feedback components. The nonlinearity is provided by 
the diodes. The voltages probed at x,y,z can be con¬ 
sidered as the degrees of freedom in our experiment. At 
fx, fy, fz external signals can be fed into the system for 
control purpose. Typical frequencies of the circuit are 
about 600kHz. 



FIG. 3. Experimental setup of the nonlinear electronic cir¬ 
cuit without the time-delayed feedback device. Experiments 
have been performed at A = llOH. 
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FIG. 2. Dependence of the period 0 on the control ampli¬ 
tude for various delay times r, from bottom to top t — T = 
— 1.0, —0.5, 0.0, 0.5, 1.0, 1.5: numerical simulations (sym¬ 
bols) and analytical expression (|^ with k = —0.9 (lines). For 
a few data points outside the stabilised regime, the period 
was estimated from a dominant peak in the power spectrum. 

Finally we have checked, whether eq.(||) successfully 
predicts the period of the unstable periodic orbit £ when¬ 
ever a few data points are accessible. To this end we eval¬ 
uated the AT-dependence of the power spectrum of the 
control signal within a regime where a periodic signal 
can be observed. Starting from r = 14.0, which differs 
tremendously from the true period, we evaluate 0 for 
K = 0.8, 0.9, and 1.0 to obtain k = —0.8 ± 0.01 and 
T = 11.745 ± 0.015 from eq.(||). The accuracy of T is in 
fact of the order of the numerical resolution of the power 
spectra. In that sense the result is striking. 


Without control the system undergoes a period¬ 
doubling cascade to chaos on variation of the resistance 
R, ending up in a Rossler-type attractor. Topological 
analysis of this three-dimensional system yielded a 
frequency of tt/T in the Floquet exponent for the unsta¬ 
ble period-one orbit of the chaotic attractor. This cor¬ 
responds to a complete flip of the neighbourhood of this 
orbit. Therefore the orbit is accessible to time-delayed 
feedback control. 

The control device consists of a cascade of electronic 
delay lines with a limiting frequency of about 3MHz and 
several operational amplifiers acting as preamplifier, sub¬ 
tractor, or inverter. The device allows to apply a control 
force of the form F{t) = —K[U{t) — U{t — r)] with r- 
range 10ns... 21/xs. Our feedback scheme consisted of 
coupling the voltage at z via the control device to fz- 

To check the coincidence with our analytical results we 
looked for periodic behaviour of our nonlinear circuit by 
sweeping the control amplitude K at hxed r. By increas¬ 
ing K the system undergoes an inverse period-doubling 
cascade ending up in a period-one state. This periodic 
state yields the desired value 0. A further increase of K 
results in a Hopf bifurcation destroying the stability of 
the periodic state (cf. fig.|l). 
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FIG. 4. Power spectrum of x{t) for t — T = 2.5ns: (a) 
below the lower, (b) between both, and (c) above the upper 
stability threshold. 

Figure ^(a) shows the main frequency at 605.5kHz and 
its subharmonic at 302.7kHz corresponding to the flip 
bifurcation at the lower stability threshold. At the up¬ 
per stability threshold a Hopf bifurcation yields an in¬ 
commensurate frequency component at 201.4kHz. We 
checked that frequency locking did not occur. By mea¬ 
suring 0 for various r and K values one obtains the 
greyshaded surface displayed in figj|. Note that for cor¬ 
rect delay time t = T one automatically gets Q = t = T. 
Within an experimental error of Ins the intersection with 
the surface 0 = r yields a straight line with r = 1.656yLts. 



FIG. 5. Measured values of the period 0 depending on r 
and K. The white plane corresponds to 0 = r. The circles 
indicate the points where spectra shown in fig^ have been 
obtained. 

Since the curvature of the surface in direction of r is 


negligible for r values close to the real period, i. e. ±10%, 
the coincidence with our analytical expression (i is quite 
reasonable in this region. Calculation of the system pa¬ 
rameter yields k = —0.31 ±0.01. For larger delay mis¬ 
match eq.(||) can still be used iteratively in the sense of 
a Newton method for detecting the exact period T. 

Conclusion - We have shown that the period of true 
periodic unstable orbits can be obtained from the prop¬ 
erties of the control signal, at least if a periodic signal 
can be realised. Our approach is based on the fact that 
the true periodic orbit of the uncontrolled system is de¬ 
formed into a fictitious periodic orbit by the control if 
the delay time differs from the true period. Our analyti¬ 
cal expression (||) relates the fictitious period 0 with the 
true period T, the delay r, and the control amplitude 
K. Peculiarities of the system enter only through a sin¬ 
gle parameter k. Of course, our result does not guarantee 
that the orbit becomes stable if the delay time is adapted 
without changing the control amplitude (cf. fig||) . How¬ 
ever, in order to keep the fictitious orbit stable during 
such an adaption process one may for example monitor 
the power spectrum of the control signal (cf. hg.^), since 
an instability is indicated by the occurrence of additional 
peaks in the spectrum. 
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